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AirrotuTxoN or on-boamd rLicurrAn HMiACBaMT 
Haim Iribargar* 

Ams Rasaas'ch Caatart NAIA> Moffacc Flald* California 


Abatract 

Tha aeacua of eoncapta and caehniqua for eba 
daaign of on>board fllghcpach unagoaaAt ayataaa ia 
cavlawad. Such ayatana ara daalgnad Co Incraaaa 
flight afficiancy and aafaty by autonating tha optl> 
■iiaclon of flight procaduraa on-board aircraft. 
Aftar a brlaf raviaw of tha origina and funetiona of 
auch ayatana, tira papar daocribaa tvo conplamntary 
Mthoda for attcdclng tha key daaign problaat 
nanaly. tha aynchaala of officiant trajoctoriaa. 

Ona Mthod optlaizaa an couta, tha othar optinisaa 
taninal araa flight; both mthoda ara rootod in 
optlul control theory. Sinuiation and flight-taat 
raaulta ara raviawad to illuatrata tha pocantlal of 
thaaa ayatau for fuel and coat aavinga. 


Par tial Liat of Svbola 

C.»C. " unit coat of fual and tim. $/lbi $/aoCi 

* ® raapactivaly 

0 ■ drag forea. lb 
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I c 
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L ■ lift force, lb 

P - dlfforantlal coat, $ par aac 

Sp^ " chruat apacific fual conauaption. par lb 

T " thruae force, lb 

t.t. - elaa. total flight, ciaa. raapactivaly. aac 
t 

t .t. - tlaa at and of cliab and baginaing of 

° ** daacant. raapactivaly 


*Raaaarch Sciantiat. 

Thk paper k dadaiod a work of tha U.S. Govcfnmant and 
thoiafore k fai the paUk domain. 
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Introduction 

Thaodore von Kdmdn la renowned for hia 
raaaarch and laadarahip in helping to aatabllah tha 
aeianclfic foundation of aaronaucical anglnaarlng. 

Aa a raaaarchar at NASA I have adalrad hia brilliant 
eoncributiona on nany occaaiona. Tharafora. I aa 
deeply honored and privilagad to ba able to pay 
hoaaga to ehla great aclantiat by dalivaring the 
lOch von Kdradn Maaorlal Lecture. 

During hia long and brilliant acianciflc career 
von Kdmdn not only eoncributad aalnancly to the 
variotia diaciplinaa of aeronautical anglnaarlng; ha 
alao played a najor role in founding aavaral of 
Chan. Thua. I would like to balicva that ware ha 
hare today ha would find aonathing of Intaraat in 
tha relatively now topic of *y lactura. which con- 
binaa alananta of parfomanca analyoia. guidance 
and control theory, and ayataa acianca. 

Pravioua von Kdmdn locturaa often praaancad a 
broad aurvay of tha lactura topic. Howaver. ny 
topic la of coo recant origin to aaka chla approach 
worthwhile, imeaad. I beliava chat the reader ia 
aarvad boat by focuaing on a few critical raaulta 
Chat ara rapraaantaclva of the currant atate of 
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knowl«ilt>2 in eha young and avolvlng flald of auto- 
laaclc fUghcpach nutnagasanc . Furchanore. I shall 
aaphaalza rasaacch conduccad primarily at Amaa 
Kasaarch Cancar during cha laat aavaral yaaca. 

Tha autoMCion of on-board fllghcpach managa- 
■anc markj cha baginning of a naw phaaa In cha avo- 
lucion of aucoaaclc fllghc concrol. For cha fine 
Cina, on-board conpucar ayacaaa will augaane or 
poaalbly avan raplaca cha plloc in planning and 
exacuClng complax fllghcpacha. This dagrae of 
on-board aucomaclon axcaada chac of axiaelng 
aucopiloc/navigaclon ayacama. which provida auco- 
oaclc guldanca only along piloc-apacifiad flighc- 
pacha. Tha hlghar lavel of aucomaclon will banafic 
boeh cha aircrafe oparacor aa Wiill an cha air 
craffic ayacom chrough incraaaod aafaey and funl 
afficiancy, and reduced plloc workload. The van- 
guard of euch on-board ayaceoa. alao raf arced co aa 
fllghc and parfocnanca aanagamanc coapucar ayacaaa 
(FPMC'a). will aoon ancar coomarclal aarvica In 
sevaral cypaa of Jec cranaporc aircrafe. 

Aucoaacad fllghcpach nanagaaanC la hero broadly 
dafinad aa compucar logic for ganaraclng a aafa, 
coaforeabla, and economical crajeccoryi on-board In 
real cloe. Thla paper praaenca fllghcpach nanaga- 
nanc cechnlquaa and algoriChma davalopad primarily 
for cranaporc aircrafe. Moraovar. cha paper empha- 
slzaa cachnlquaa chaC have been evaluaCed In pllocad 
almulaclon and fllghc caaCa and are being Impla- 
menced in commarclal ayacama. 

Becauaa Incareac in aucomacad fllghcpach man- 
agamanc haa bean moclvacad primarily by incraaaing 
fuel coaca. cha focua of raaearch haa been on find- 
ing compuCar-lmplamanced aoluclona Co eha minimum- 
fuel and coac-crajeccory problema. Thsrafore. Che 
main purpose of Chla paper la Co daacrlba currencly 
uaad algorlchaa for on-board calculaclon of fuel and 
coat efficianc crajeccorlea. Also, Cha Incarface of 
cv>^ algorlchm wicb plloc diaplaya and ocher guidance 
ay. cams will be reviewed, wlch raferanca co an 
Implamcncacion racencly avaluacad in fllghc caeca ac 
Ames Raaearch Cancer. Finally, rasulcs from almula- 
clon and fllghc ceacs will llluacraee Che efficacy 
of chase syscema co opclmlze cyplcal airline fllghc 
mlsalona . 

The paper begins wlch an ovarvlaw of cha func- 
clons and general scruccura of fllghcpach managaoMnc 
ayaCema. To simplify on-board ImplamancaCion, Cha 
analysis is divided Inco cwo complcmancary problems: 
en rouce and caivlnal area fllghc. Each of cha cwo 
problema la examined In saparace aecclona, which arc 
complace In chamsalvas In chac chay Include cha 
derivaclon of che on-board algorlchm and almulaclon 
or flighc-casc reaulcs. Tha sacciona are organized 
CO ssclsfy cwo classea of readers: Chose Inceresced 
In Che analysis and on-board ImplamancaCion who will 
find analyclcal decails aufflclenc co cranslaCa Cha 
algoriChma inCo compucar coda; and chose seeking a 
quick overview, who can skip che analyala and con- 
cencrace on che Incroduccory and reaulcs subsacclons 


Fllehcoach Manaeeiaenc FuncClona and Problems 

A convenclonal aucoplloc la dealgnad co crack 
various cypaa of fllghcpacha, Che slmpleac of which 
are holding a specified alclcude, speed, and head- 
ing. More complax casks performed Include cracking 
of chraa- or four-dlmanalonal curved crajeccorlea 
(che fourch dlmanalon la clme). Tha almpllfied block 


diagram of a fllghcpach Cracking aucoplloc Is 
ombadded In Fig. 1, Ics kay elesanca are a compan- 
sacor module, a aansor/eaclmacor/navlgaclon syseam 
nodule, and a summing Juncclon. Errors beewaan 
cooaanded and aircrafe scacas are conclnuously 
nulled by cha acclon of cha feedback loop. Hareco- 
fora, command Inpuca co che aucoplloc have bean 
genaracad by cha plloC, buc currant devalopmancs 
ara changing this process. 


KEAL TIME WPUTt: 

• MNFORMANCE CniTEgKM 

• wiNoa 



Fig. I scruccura of fllghcpach managenanc syscom. 

Aucomaclon of fllghcpach Dunagomanc la cha 
process of ganeraclng InCalllganc command InpuCs Co 
cha aucoplloc by an on-board compucar, as lllus- 
craced In Fig. 1. Tha data base for such a syseam 
la axCenslva. IC Includes a decallad model, often 
scored In mule Idlnans Iona 1 Cables, of aerodynamic 
and propulsion syscem performance. Ocher aleowncs 
contained In Ic ara alrllna roucaa, terrain and 
alrapaca models for carmlnal areas, cha aircraft 
operational envelope, and fllghc rules. In auammry, 
cha compucar muse have available che same kind of 
informaclon required by a plloc for safe and affl- 
clenc aircraft operaclon. The icema llacad under 
real elae Inputs In Fig. 1 ara so dafinad because 
chay are frequently updated during flight alchar by 
che pilot, che navigation system, or by a data link 
CO a ground facility, such ao an air Craffic control 
center. Perhaps ic Is surprising chat che perfor- 
mance objective la included as a real-clma Input; 
however, many conditions can arise that will require 
It CO be changed during fllghc. 

The decs base and Che real clma Inpuca ara 
operaced upon by Che algoriChma In Cha fllghcpach 
managaieenC compucar to generaca efficient and 
conflict-free crajeccorlea. Thla process la analo- 
gous CO cha work of che flight crow. In chac it 
Involvaa planning, monicoring, and revising the 
crajactory throughout che fllghc. Traiaccorlea are 
first synchaalzad (planned) In "fasc elms'' 
chac la. In a elms Interval chac la a small fraction 
of che actual fllghc clme. In this crucial scap, 
cha algorlchm computes che entire future fllghc 
history from che current position co cha landing 
point. Than Che computed ccajeccory la aCorad and 
finally transformed Inco real Clma cramand Inputs 
for Cha aucoplloc. During flight the syaCsm monl- 
Cors cha trajectory for Inclplanc confllcCa with 
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Intruding aircraft and for axcaaalva tracking arrora 
cauaad by uiwodalad wlnda and othar diaturbancaa. 

It alao nonitora pilot inputSt auch aa changaa in 
tha parfocuanca crltarion or daatlnation point. 

Thaaa and othar condltiona can trigger a rovialon 
or coaplata recalculation of tha trajectory by tha 
faat-tlM algorltha. Laat anyona trlaa to acquira 
a ayataa with thaaa capabilitiaat I haatan to add 
that aoM of tha "aBart" functiona anvlalonad hara 
are not yac available in tha currant ganaration 
flightpath nanagaMne conputara. Howavar, raaaareh 
la rapidly noving the atata of the art coward chair 
raaliaaCion. 

Of tha functiona outlined above, on-board 
opciBlaatlon of crajaccoriaa haa racaived tha aoat 
attention, since it liaa at Che haart of the flight- 
path nanagenanc problaa. In airlina oparationa it 
ia generally agreed that tha aoat uaaful parfomanca 
criterion la tha total coat, J. of a aiasion, which 
ia defined aa tha aua of fuel coat and tine coat, 

J ■ CfWf -f C(tf, with flight ciaa tf unapaclflad. 
Mlniaua fuel and aininua tiaa criteria are apaclal 
caaas obtained by aatcing C( or C{ to zero, 
raapacCivaly. In ordar Co aaat an ai.'line flight 
schedule or an aaaignad landing ciaa alot, ainiaia- 
Ing J with apaciflad arrival ciaa is alao of 
intaraac. 

It is convaniant to separata trajectory prob- 
laaa into two clasoas; naaaly, on routs problaas 
with fllghtpaths longer than approxiaacaly SO n. al. , 
and caralnal area probloaa with paths shorter chan 
SO n. al. In an route flight, the patha are pra- 
doainancly long sections of straight lines with a 
nagllgibla percentage of the flight tlai Spent in 
turns. Thus, turning dynaaics can be naglactad and 
only vertical plana dynaaics need to be aodalad in 
opclalzing tha an routa case. This problea ia 
studied in the next section. In caralnal-arsa 
flight, vertical and turning aanauvara tend to occur 
siaultanaoualy and in comparable ciaa intervale. 

Thua the dynaaics of both types of aocion auat be 
aodalad in trajectory opclalsation. This aore dif- 
ficult probloa is studios last. Although solution 
to both problaas have been carried into slaulation 
and flight tsacs, tha aorging of the solutions 
required in a full-alsslon flightpath aanagaaanc 
syscoa roaaina to be sccoapllahed . 


where is tha coaponant of horlsontal wind 
valocity along tha ground track diraetion. Tha 
quantity V,, can ba a function of sltituda, but 
dyntalc offsets of wind shear aa wall as tha verti- 
cal coaponant of the wind do not play a significant 
rols hara and ara naglactad. In airplanes, unlike 
in aoat types of aisailes, aass flow owing to fuel 
bum is relatively alow and does not naad to ba 
aodalad by a atata aquation. Instead, tha slowly 
changinit mss of tha aircraft will ba traatod aa a 
tiaa-varying paraaetar. 

In this paper th» aquations of notion ara 
furChar ainplified by co^inlng altitude and air- 
apaad into a olaglo atata variable, specific energy: 

E - h ♦ (l/2g)V* (5) 

Diffarantiating Cq. (S) with raspact to ciaa and 
subacltuclng Eqa. (1) and ( 2 ) into Eq. (5) yields 

dE/dt > (T - D)V/W ? k <6) 

The control variables in Eq. (6) are airapaod V 
and thrust T, or its related quantity, throttle 
sotting V. This ia the ao-callad anargy-stata 
oodal, which has bean widely used in trajectory 
optiaiaation problaas.^"’ Its utility dopanda 
entirely on tha naturo of the application. For tha 
quasi-scaady crajaccoriaa coaaonly found in clisd>, 
crulaa, and descent of transport aircraft, the 
anorgy-etaca aodal provides especially aiapla 
on-board algorlchas, an wa shall saa. 

Optiaal Control Foraulation 

In a previous section, the aoat iaporcant por- 
focaanca criteria that arise on-board flightpath 
aanagaaanc ware anuaaraeed. These criteria will be 
shown hara to be assancially equivalent when forau- 
lacad aa probloaa in optiaal control. 

Consider first the ainiaua-cosc criterion with 
a specified range to fly and no explicit constraint 
on flight tisw. This criterion can be written as 
an integral coat function: 

r'f /»^f 

J - J (CjWj + C^)dt = J P dt (7) 


En Route Fliihtpath Optiaiaation 

The point naas equation of notion for flight in 
r.lie vertical plane can be written as 

dV/dt - g(T - D)/W - g sin Y (1) 

dh/dc • V sin y (2) 

dx/dt ■ V cos Y (3) 


With Eqa. (4) and (6) as Che state equationa, the 
Hoailconian of optimal concrol** is 

H - Cjfilj ♦ Cj + i|/g(T - D) J + ii^(V + V^) (8) 

where and are the coatatea. On an extrsaua 
crajaecory the Hanilconian achlevaa its nlniaua 
with respect to the controls V and T, and the 
costates obey the linear dlfferancial aquations 


In normal flight aanauvara of transport aircraft the 
flightpath angle rates ars such chat yV/g <c I and 
|y| i 10*. These conditions allow us to take flight- 
path angle as a control variable with lift calcu- 
lated fron Che constraint L - W cos y and to 
approxlasce the cos y factors in Che above equa- 
tions by unity. The effect of a horizontal wind, 
when its aagnltude is a saall fraction of the air- 
speed (one third or less), can be included by aodi- 
fying Eq. (3) as follows: 

dx/dt - V + s (4) 


HCfVf + 
3E 






or 


^x-° 

- constant 


( 10 ) 


Since specific energy and range at the final cine 
are specified, the values of the coatatea at Cf 
arc free constants, which are used to deteraine the 
desired final states E; and d£. Tha proper choice 
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of Ch*B« conatMCa coapriaaa cha waLl-knuwn cwo- 
polnc boundary valua problaa of opriul concrol.^ 

Tha aolucion of tlita problaa will ba addroHaad in 
tha naxc aactiont Bacauaa cha final cltaa ia uncon» 
acrainad, cha Hanilconian alao obaya cha cundicio'.', 

H - 0 (11) 

for all c. Noca chac <|<£ ia a funcclon of cIm 
wharaaa ia a conacanc Indapandanc of cImi 

Howavart it ia noc nacaaaary to inCagrate Bq. (9) 
along the Crajactory Co obCain iig(E). Equation (11) 
providaa a firaC ineagral and whan coabinad with 
Eq. (8) yialda aigabraically, 

Conaidar naxc cha aaaa parfominca crltarion 
aa Eq. (7) but with cha final CiM c^ apeclfled. 
Thia changaa only Eq. (11), which bucoaaa 

H - (12) 

wharo ia a conacanc for aach crajactory. But 

for each conacanc c)^ obtainad for u particular 
fixed-cima opClaMS crajactory, ona can dafina 
C( - C(. - Cc and conaidar tha aaaa trajectory nlao 
CO be a aolucion of Cha fraa-ciaa problan, with 
replaced by C{.. Obvioualy, the cima-coac factor 
cannot be apaclfiad indapendencly in Che flxad-clM 
problan. Thua, the fixad-cine problan wich apaclfiad 
range Mkaa phyalcal aanae only for Cha ninimn fuel 
perfomanctt index, ■ 0. Moreovar, every fixed- 
cinie, nlnlimim fuel problan can be fomulacad aa a 
free-cine problan with the clne-coac factor 
choHen as co achieve cha desired final time. 

Finally, consider Che ninlnun fuel perfomance 
index, Cc ■ 0, with fixed final tine and no con- 
straint on range. This problem occurs in generating 
mininun fuel delay maneuvers. Here the cransvaraai- 
ity conditions of optimal control require i’ 0 
(Ref. 4) and a particular solution will generally 
yield a nontaro tl, siiy C(, and some range 3f. 
However, this solution will be identical to the^ 
f roe-clsHi, _f ixed-range problem in which ■ -Cj 
and df • dj. 

The esaencial equivalence of the various prob- 
lens implies we can concentriice, without loss of 
generality, on slgoclchns chut solve the free-final- 
time problem. 

Solution Baaed on a Simplified Approach 




+ j’^‘^(P/|Ell^^,)dE (14) 

®f 

whara E^ and Ef are Che glvan Initial Uilnb and 
final descant energlas, raspactlvalv. Ths crans- 
fomaclon uses the asaunpclon chac the energy 
cliangas aonotonlcally In cha climb and descwnc. 

This places strict inequality conscrainta on E, aa 
shown in Eq. (14). Alao in Eq. (14), the incagra- 
tion limits hava been reversed in the dascent cost 
term. In this formulation tha coaC function la of 
mixad form, containing two integral cost taima and 
a tacminal coat term contributed by the cruisa 
segswnt. 

Wich cha change in indapandanc variable from 
time to onergy, the state aquation for apacific 
energy is eliminated, leaving Eq. (4) aa tha only 
state equation. FurCharmora, we noca chat tha 
performance funcclon (Eq. (14)) depends on the dis- 
tance state X only through the aia of cha final 
vaiuea of climb and daacent distancaa d^. ■*- d^^^. 
Therefore, the state equation tor Che dlatancc is 
rewritten in terms of this stus: 


We now introduce the approach of Ref. 3 by 
assuming chat the trajectories are composed of three 
segments; namely, .'i climb, a cruise at constant 
specific energy, and a descent, as illuscrated in 
Fig. 2. The cost function (7) can then be written 
as Che sius of che costs of the three segments. 



climb cost cruise cost descent coet 


where \ designstes che cost of cruising sc a given 
specific energy E^. Next, we transform the inte- 
gral cost terms in Eq. (13) by changing che indepen- 
dent variable from time co aoeclfic energy, using 
che transformation dc • dE/E: 


♦ <*4. * 'i.. 


U5I 


Here, che transformation dc “ dE/E was used again. 
Also, Eq. (IS) provides for Independence in che 
apaciflcaclon of che climb and dascent wind velocl- 
clee Vyup and Vyj,^. Generally, different wind 
condlclone will prevail In phyalcaily different 
locacione of climb and descent. The wind velociclee 
can also be altitude dependent. The effect of 
altitude-dependent winds on cha optimum trajeccorlea 
is dlacusaed in Ref. S. 


Necessary conditions for cha minimisation of 
Eq. (14), subject to tha state equation (Eq. (15)) 
yield cha following exprasalons for che Hamiltonian 
and coacaca equaclona, raapectlvely: 
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dip/dE - -(3H/a(x^p + - 0 


( 16 ) 

<» 7 ) 


The HakllConian la aininlaad with raapact Co 
two palra of control varlablaSi ona pair appllcabla 
Co cllab (Vup and ir„-) and cha othar co daaeanC 
(Vj„ and Noca'^chac chroccla-aaccing it 

rachar chan chruac la uaad aa a control varlabla. 

In genaraL. chruac and fual flow ara nonllnaar 
Cuneciona of <r aa wall aa of alcicuda, Mach nun- 
bar. and caaparacura. Slnca each cam undar cha 
■Inlnlzacion oparacor In Eq. (16) concalna only ona 
of cha cwo palra of concrol variablaa, cha ■Ininiaa- 
clon alnpllfiaa Into Cwo indapandant mlnialsatlona , 
ona Involving cllnb concrola. cha othar. daacant 
ooncrola. Alao. slnca Cha rlghc-hand aide of Cha 
coacace aquaClon (Eq. (17)) la zaro. la conacant. 


Naxc wa axanlne ch« CranaveraallCy condiClona 
applicable Co this fortnilaclon. The baalc con- 
scrainc In chls problem la chac che range of cha 
crajeccory be df. Howaver. df la a paramacar In 
Che cranaforud cose funcclon. Eq. (14). and noc a 
scaca V:i;‘labla. The final value of Che scacc 
varlabla dup -f djn In chla focbulaclon, sub- 
Jacc only co cha Inaquallcy conacralnt d^p-f ddnldf ■ 
Thla conatralnc la, of course, naceaaary for a 
phyalcally meaningful raaulc. The inaquallcy con- 
acralnt can be handled by solving two optimization 
problama. one compleCely free (dyn * dj„ ^ df). the 
ocher conaeralhad (d^p -f d<in ■ df). and Chen choos- 
ing cha crajeccory with cha lowest cose. Physically. 
Che comparison Is between a crajaccory wlch a cruise 
segment, and one without a crulsa segment. Consid- 
ering first Che f rea-teralnal-acata case 
dup -«■ dju < df. we obtain the following relation for 
the final value of che coatace 


9(d, - X _ - X j_)X 

) . i SIE ±L- 

3(x X. ) 

' Up dn' 


®"^'’‘up"“up'*dn’^dn 


(18) 


In both cases i;ha choice of a coacace dacar- 
mlnas a particular range. Since tha functional 
ralatlonahlp batwaan Chase variables cannot be 
decaralnad In closed forn. it la naeassary co itar- 
ace on cha coatace value In order co achieve a 
specified range df. 

Tha laac nacasaary condition appllcabla to 
chls formulation Is obtained by making use of cha 
fact Chat Cha final value of the time-llke Indepen- 
dent variable E Is free. Its final valua la cha 
upper llmlC o' Incagraclon Ec In Eq. (14). Appli- 
cation of the fraa-flnal-cima CranaveraallCy condi- 
tion In Ref. 4 provldas cha following condition: 

(h + OUdf - d^p - dd„)^<E)I/9E))B.£ - 0 

( 19 ) 

which, whan evaluated becomes 

IH + (dc(dX/dE)))j._g^ • 0 (20) 


whara d^ Is che crulsa discanca. 

Pa)uaclon (20), cogachar wlch knowladga of che 
sallanc charactarlstlca of cha crulaa coat > and 
cha Hamiltonian H, can be used to dataralna the 
atructural dependence of cha optimum Crajactorlas 
on che range. 

Crulaa coat at a crulsa energy E^ and crulsa 
spaed Vq Is computed from Che relation 

X(E^.V^) - P(T,E^.V^)/(V^ + V^) ; 


constraints: < 


where tha denominator Is cha ground spaed In the 
fllghcpach direction. Examination of cha term 
containing X In chu relation for che performance 
function (14) shows chat the valua for X should 
be as small as possible at each cruise energy In 
order to minimize the total cost J. Tharafore, Cha 
crulsa-spaed-dependance of X la allmlnacad by 
minimizing the right side of Eq. (21) with respect 
to V^: 

X(E^) - min P(T,E^,V^)/(V^ ♦ V^) (22) 

'^c 


This la the cransvert.:llcy condition for cha fraa- 
fin&I-acace problem wlch catmlnal coat.** It shows 
chat che constant coscata value la Che nagaclva of 
che cruise cost. 

Next, consider Che case of trajectories with 
no cruise segment. Tlten. the middle term of 
Eq. (14) drops out and the performance funcclon 
concalna only the Integral coat terms. Thla is Che 
case of che specified final state df ■ dup djq; 
the corresponding CranaveraallCy condition yields 
<ji(Ec) “ constant. In practice It is not necessary 
to compute tha constrained cerminal-stace crajaccory 
if a valid fraa-carminal-scaca crajaccory exists, 
that Is, one for which df > d„p dj„, since tha 
addition of a terminal constraint can only Increase 
the cost of Che crajeccory. Therefore, this case la 
noc considered further In this papar. 


In chls papar, X and V,. arc alwayc. aasumod to 
be cha optimum cruise cost and cruise spaed, respac- 
clvely, sc a particular cruise energy E^. 

Except In high wind shear, tha cruise cost as 
a function of crulaa energy exhiblta cha parabollc- 
llke shape shown In Fig. 3. For subaonlc transport 
aircraft, che minimum of Che crulsa coat wlch 
respect co energy occurs cloae co the maximum energy 
boundary. This characteristic of cha cruise cost 
prevails for esaentlally all values of cha perfor- 
mance function parameters Cf and C(. The quanti- 
ties defining the opclmta cruise conditions are 
Ecopt •»“* ^opf Eq. (20), the derivative of che 
crulae-cosc function multiplies the crulsa discanca. 
Except under axersme wind-shear conditions, che 
derivative is monoconlc and crosses the zaro axis at 
Ec - Ecopc By observing Chat climb and descent 
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Fl|. 3 Crula* C9tC Function. 


controls do not occur slaulcsnaously In tha taraa of 
Eq. 06)f and by subatlCuting Eq. (18) In Bq. (16)i 
we can aaparata H Into climb and dascant compo- 
nants as follows: 


whare 




'up • 

up 

%p 


p - X^up * 

L *'fc. J 


I ■ min 

'^dnL I^I1 b<o J 


'dn 


(23) 


<24) 


Nuaarical studlas of Eq. (23) for savaral 
Bodals of subsonic turbofan aircraft show 
H[Eci<\(Ec)) > 0, for < Ecg-^. Consider first 
tha casa In which H(EciX(Eq)] > 0; than Eq. (20) 
can ba solved for the cruisa distance d^: 

d^ - -H(E^.A(E^)l/(dl/dE)g,j.^ (25) 




Fi|. 5 Energy vs range. H ■ 0 at B^. 
Thrust Optimiaation for Mtnlaai Fusl Trsjsctorias 

Evslustlun of tha Hsadltonlan would bs simpli- 
fied if one of the two pairs of control vapisblas. 
airspaad or throtcla. could ba aliminatsd a priori 
from tha minimisation. Slaea the pair of throttle 
settings ?up *dn thought to ba naar its 
limit, wo shall look for conditions whara extrema 
settings of the throttla are optimum. Wo axamins 
here only the miniaum fuol case C{ ■ 1 and 
C( ■ 0. with winds sot to saro in order to simplify 
tha derivation. However, tha results can ba 
extended to tha casa in which C( d 0. 

For minimum fuel performance, tha two tame in 
the Hamiltonian Eq. (24) become 


Since dl/dE < 0 for Eg < Ecopt* d\/dE 

approaches zero as E^ approaches Ecoac. the 
cruise distance must increase without limit as Eg 
approaches E^g.^. Although numerical calculations 
show chat the value of H tends to decrease as E^ 
increases, the race of decrease in dX/dE is more 
rapid and doninaccs the behavior of d^. Figure 4 
shows the family of crajeccories. obtained if H > 0 
for all values of Ef In this casa, Incerastingly, 
nonzero cruise segments occur at short ranges and at 
energies below the optimum cruise energy 

Consider next Che case in which 
H[E(.,X(Ec)] ■ 0" Then, d^ ■ 0; that is, no cn'-'-ie 
aegmanc is present for dX/dE < 0. However. 

Eq. (20) allows d^ to ba nonzero if dX/dE ■ 0. 
This implies that for H ■ 0, cruise flight is opci- 
Bum only at the optimum crulae energy Ecopt* 

Figure S shows the family of crajeccories tor this 
casa. 


1 


up 


nin 

”up''^up 


up 


*^dn 


min (26a) 

”dn’''dn 



An accurate model for thrust and fuel flow 
gensrally includes the functional dapendeneies , 
T(«.V,h) and U{(v.V.h). In addition. Chose func- 
tions must be corrected for nonstandard camperaturas 
and bleed losses. 
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In •arllar work on aircraft trajactory opclBi- 
iation«^ a alaplar aodal for fual flow and chruat 
waa of tan uaads 

li, - TSp^.(V,h) } T^in(V.h) S T i T^^<V,h) (27) 


Tha critical aaauaptlon in Eq. (27) la indapan- 
danca of tha apaclflc fual conauaiptlon Spc froa 
thruat. Tha vlrtua of thla nodal llaa In tha 
Inalght it ylalda into tha alnlnua fual problan. If 
Eq. (27) ia aubacltutad into Eq. (26b) • ona obtalna 


up 


»FC“ - <»=FC>'' J 

V I Jt 

up 


’ '^dn '^dn - J 


(28) 




For any flxad valuaa of Vup or Vdn> tha 
operand functlona Kup and K^n with ratpacc to 
chruat ara hyparbolaa with polaa at T„p - D and 
Tdn " D, raapacClvaly . Tha nunaracor zaro with 
raapacc to chruat nuat ba to tha laft of tha pola on 
tha chruat axla for apaclflc anarglaa laaa than 
crulaa anargy. Thla inpllaa that nax^num Chruat 
nlnlalcea K^p and Idla chruac alnlnlzaa for 

any B < Ect and provaa chat tha Uniting valuaa of 
chruac ara optlnun for chia propulaion nodal 
throughout tha cllnb and dbaeant trajactorlaa. Thla 
raaulc alao Inpllaa chac tha daparcura fron cha 
Uniting chruat valuaa found for cha nora ganaral 
propulaion aiodel la directly accrlbucabla to cha 
nonlinear dapandanca of fual flow on thruat. Con- 
varaaly, the need for throttla-aactlng optlnlzaclon 
can be datamlnad a priori from cha fuel flow veraua 
thruat dapandanca for a particular engine. Such 
data are found in the engine nanufacturer’a parfor- 
nance handbook.* 


Effect of Model Characterlatlca 


We have aaan In a preceding aactlon chat the 
value of cha Hanlltonlan conputed at crulaa energy 
Ec dacamlnea the acructure of cha crajeccorlaa 
naar crulaa. It la poaalble to relate the exlacenca 
of crulaa below E^gpe to apaclflc angina and aaro- 
dynanlc nodal paranacara. This la dona by aubatl- 
cuting truncated Taylor sarlaa axpanalona of fual 
flow and drag aa functlona of alrapaad and chruac 
Into cha axpreaalon for the Haallconlat.. Tha loca- 
tion of Che alnlnua with raapacc to tha controla aa 
well aa cha value of H can Chan ba dataralnad aa 
functlona of Cha Taylor sarlaa coefficients at 
£■£(.. Theaa and other calculaClona are carried 
out in Ref. 3, Hare the reaulcs ara suant.rlzed for 
two types of angina characteristics. 

Casa Ai Sfc Independent of Thruat . Tha 
structure of cha trajactorlaa for cha case in which 
specific fuel consunpclon is Indapandant of thrust 
Is given by cha fanily of trajactorlaa In Fig. 5. 

In these trajactorlaa no crulaa sagnant occurs 
unless cha range excaada a cartaln nlninun range, 
and chan crulaa takes place at cha opclnun crulaa 
energy Eggp^. Furchamora, It can ba shown that 
cha opclnun cllnb and descant spaeda ara equal to 
the crulaa speed at the crulaa energy.’ Thus thara 
is continuity of cha opclnun speeds at entry and 
exit of the crulaa sagnent. If It la present, or at 


tha cranaltion fron cllnb to descant, if it la not 
prasant. Tha continuity of the apaada slnpllfias 
tha design of cha Interface beewaan thla algorithn 
and tha guidance systan for flying the trajactory. 

Fravloualy It waa ahown that for cha chruac- 
Indapandant caaa, tha opclnun cllnb thruat la 

tha naxlnun thrust, and cha opclnun daacant chruat 
la cha nlninun or idla thrust. Thaaa and cha above 
alnpllfylng characteristics enphaslse cha Inporcanca 
of cha asaunptlon underlying Casa A. It auggasca 
chat even If cha aasunpeion la not conplaceW satis- 
fied, ona should avaluata Its efface on the opcl- 
nallty of the crajeccorlaa. 

Casa B; Sar Thruat Dapandant . Figure 6 plots 
Spc and fuel flow as a function of thrust for a 
typical In-sarvlce turbofan.* Over the conplata 
thruat range, S|p(. Is aaan to be acrongly thrust 
dependant, approaching infinity at low thrust 
values. However, at typical c.lnb or cruise 
chrusca, corresponding to cha upper half of tha 
chruac range, tha variation In Spc is only about 
5X. Tha dashed line through the origin glvas cha 
bast constant Spc approxliaaclon to cha fuel-flow 
function. Conparlaon indlcataa an excallant natch 
at high thruat, but significant errors at low 
thrust. For sons applications, tha constant Sy^ 
aasunpeion of Casa A may still ba adequate, If 
errors In fual flow at low thrust are considered 
unlnporcnat . 



Fig, 6 Syc and fual flow vs thrust; typical 
turbofan. 

By using £n engine modal slnilar to that 
defined In Fig. 6, it was shown in Raf. 3 that if 
thrust is a free but bounded control variable, cha 
opclnun cllnb and descant thrusts and airspeads 
converge coward the opclnun cruise thrust and air- 
speed as tha cllnb and dascant energies approach 
cha cruise energy. Thla result appllaa to all 
cruise energies, including those less chan tha opcl- 
muai crulaa energy Ecopc acructura of cha 

crajactorlas is sinllar to chat shown in Fig. S, 
except chac at cha maxlnun (or cruise) energy the 
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traj€CCori*a will h«v« a rounded rochar than a 
paakttd .ippaarance (bacauaa of tha continuity of 
thruat at BaxlauB anargy). Hlnlnixatlon of H by 
coaputar haa ahown that tha thruat ia aaxiana or 
Id la for about tha lower three fourtha of the energy 
range batwaan Initial (or final) and BaxlMai energy 
and t'.wn convargaa to tha cruiaa value aa cruiaa 
energy ia approached. Thia corroboratea the theory 
that the optiaua cruiaa energy ia approached aayap- 
totlcally aa the apeciflad range incraaeea without 
bound . 


K different reault in obtained if thia fuel- 
flow modal in uaad in conjunction with the conatralnt 
that thruat be net to ita maximum value in climb and 
CO ice minimum value in daucent. Such a conatralnc 
might be invoked in order to aimplify the ainlmiaa- 
cion of Che Hamiltonian, and it given rine to non- 
zero cruiaa aagmanen below the optimum cruiaa energy. 
The acruccure of Che Crajeccorlen in thia caae ia 
llluacraced In Fig. 4. Again, optimum cruiaa energy 
ia approached aaympcotically for Urge rangee. 

Computer Alaorithm . The climb and deacenC pro- 
filaa are generated by integrating the atata ei^ua- 
tion (IS) from Che initial energy to the 
maximum or cruiaa energy E^. For thia purpoaa, 

Eq. (15) ia aeparatad into ice climb and dencenC 
componenta aa followa: 


dx 


- (V^p coa Yyp + V^p)/e]e>o! ^ 

Ydn W^'^Oexu- ’‘dn(^f^-«i 


dx. 

-ia . (V 

dE ^'^dn 


H29) 


with E ■ (T - D)V/W. Flightpach anglea are not 
defined within the reduced dynamlca of the energy 
atato model and were previoualy uanumed to be amall. 
Hevercheleaa, during the Incagratlon of the trajec- 
tory, the flightpach anglea for climb and deacent, 
V^p and Ydn' computed by uaing valuen of 

altitude and dlatance from two aucceaaive energy 
levela. The uaa of theae computed flightpach anglea 
in Eq. (29) increases somewhat the accuracy of the 
climb and deacent dlatance Integraciona. In addi- 
tion, equations for time and fuel ace also inte- 
grated during climb and deacent. 


At each energy level the optimum airspeeds and 
thruat settings arc obtained aa the Values that 
minimize the two components of the Hamlltonlen in 
Eq. (24). Tht minimization of ths Hemlltonlen ia 
carried out by the Fibonacci search technique. ' 
Fibonacci saarch la basically a onc-variabic mini- 
mization procedure. It ia adapted here to two 
variables by applying the tachniqua to one variable 
at a time, while holding the other variable fixed. 
Convergence to the minimum ia achieved by cycling 
between the two variables several timaa. Prior to a 
search over a given control variable, the limlta of 
the regiona for K^p and are computed In order 

to keep the search interval aa small aa posaibls. 

As previously explained, the choice of X in 
the Hamiltonian decetmines the range of the trajec- 
tory, but the exact functional dependence between X 
and range cannot be determined explicitly for the 
various weights, wind profiles, and other parameter 
values encountered in real-time operation. Itara- 
clon on X must therefore be used to achieve a 
specified range. Since each itaration step requires 
intagraclng Eqa. (29), it ia important to mlnimlza 


the number of iterations. This ia accomplished by 
updating tha aaclaece of X at each itaration 
btepi uaing an approximation of tho functional rela- 
tionahip between X and d(, ee illuacreced in 
Flgt 7 for a typical trenepnrt aircraft. For any 
range df < dai„i the astimece le computed from 

X - K/df •¥ B (30) 

where the eonetente A end B are updated after each 
trajectory Integration, Further details of this 
procedure era given in Ref. 7. 



Fig. 7 Typical cruise coat va range relaclonahlp 
for example aircraft, 

If in addition to range tha time of flight ia 
slao apacifiad, a second iteration loop involving 
the time-cost factor must be implemented, This 

loop usee a procedure analogous to the one in the 
range itaration. 

Anothar elmnent of the algorithm compansetae 
for tha weight change ceueed by fuel bum. The 
effect of the change in weight on tha optimum tra- 
jectory is estimated by two mechode. The first 
merely updatea tha weight in the calculation of E 
during climb and deacent. This enauras Chet Che 
crejectorlee generated by intagracing Eq. (29) era 
baaed on an accurate model of the elrcrefc. 

The second method atcampea to correct the 
optimiaacion of the trajectories by estimating tha 
change in the cruise coat X. rhla ia done by 
using the weight of the aircraft at the end of 
climb, Chat la, at energy E^ co compute the value 
of X, It le imporcent to use the weight at 
rather then the weight at E^^ to compute \ 
because Che acnaitivlCy of the optimum controls to 
cbmngea in \ increases aa Che aircraft energy 
approachaa E,.. The fuel consumption for Che ontire 
climb trajectory, F^n, is eacimated at the start of 
climb from the empirical relation 

"up-",", - <’•> 

where is an aircraft-dependent conatant, and 

*^ref ia a typical Initial climb weight. This rela- 
tion aeclaecaa the climb fuel weight to about lOX 
accuracy, which la adaquece for this purpose. Siml- 
lerly, the weight at the end of cruise, if a cruise 
aagmsnt ia present, is used to compute X for the 
descent optimization. Tho cruise fuel conaumpr^on 
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Fc Is dstsnilnsd fro* tha ralaclon 

Fg - (32) 

whara W it cha avarsga (ual-flow rata and V« tha 
avaraga ground spaad during crulaa. Tha calculacion 
of W and ia daacribad in liaf< 1> 

Tha coaputar taplaaancaclon iricludaa both tha 
fraa- and constrslnad-thrust caaaa. For the 
const ralnad-ehrust caaa, tha crulaa dlstanca la 
coaputad (roa Eq. (25) • Howavari bacauaa dX/dC 
approachas xaro as approachas Ecop^i thnrti la 
a practical Halt to cha use of Eq. (25) • datar- 
nlnad by cha nuaarical accuracy of coaputlng dX/dE 
for Ec in cha nalghborhood of Eqo„(. A practical 
Liaic for Ec value for which 

\ ■ l.OlXopc. The coral ranga of Cha trajactory 
obtained for this value of \ ia rafarrad to as 
‘^aax* ^31 crajectoriae requiring longer ranges 
than dgi,^ ara aaauaad Co crulaa at E^opc and to 
contain crulaa aegaenCa of length d^ ■ df - dup - d^ui 
where dyp and dgn coaputad for X • l.OiXopc. 
In cha fraa-chrusc caaa, nuaarical difficulciaa can 
arlac In alniatcing Eq, (24) an E^ approaches 
Ecopt* value of l.OUop^ has also bean found 

to serve as a practical criterion for coaputlng Che 
longest range without a cruise sagsanc at E^gp,.. 

A computer pv-ograa of the algorltha has been 
impleaented in I'ChTRAN IV and la described in detail 
in Ref. 7, There are approxiucaly 2,400 FORTRAN 
instructions in the prograa. 


for a fixed final tiae the fuel conauapclon 

aa read froa cha ST locus aust bt less than cha 
fuel consuaptlon as read froa tha SR locus. More- 
over, at c,£ tha corresponding range dj on the 
SR locus aust be greater chan tha range dg on cha 
ST locus. Thaae charactaristica follow froa cha 
opclaality of the loci. Aa the claa-coac factor 
Ct is changed over a range of poaitlva end nagaciva 
values with ranga held fixed at di, the fual-tiaa 
locus MC(di) of all ainiauai cost, fraa-ciaa tra- 
Jactorias is generated , Tha ainiaua of this locus 
with respect to fuel consuaption is attained at 
Ct ■ 0 and contributas one point on Che SR locus. 
Foinca to the laft of cha ainiaua correspond to 
C( > 0 and points to cha right Co < 0, Froa 
arguaanCa in a prticeding section (Opciaal Control 
Formulation), the HC(d|) locus aust have one point 
in coason with the ST locua; naaaly, at tha ciar 
Cad Fig S. Furcharaora, cha two loci aust ba 
tangent at that tine and can have no other points 
in coaaon. 

Aasiaia that the specif led final tiae t) is 
such chat Caf S t; i tag. Tha dlffarance 
Cl ” Caf is the delay with respect to Che ainiaua 
fuel final tiae. Mao asauae that d^ is cha 
shortest distance froa Che current aircraft position 
to Che destination point. We now as, whether a 
screeched pach, say dj > d) cun give lower fuel 
conauapclon Chan d^. But ic follows from Che rala- 
clonship illustrated in Fig. 8 chat NC(d^) HC(dj) 
for d} > d( and Cj £ C^« Thus pach stretching 
is noc optlnuB for t^ if Cg|j, 


Optiauin Tiae-Delay Trajectories 

Many air travelers have experlanced in-flight 
delays ac cheir dascinacion airports. Such occur- 
rences lead CO the iaportant flighcpach managaaanc 
problem of ninlaizing cha fuel loss for u specified 
delay in the landing ciaa. Assuaing ic is not fea- 
sible to absorb the delay on Che ground or ac an 
ultemaCive airport, we can divide delay manauvera 
into two cypea: slow-down and paCh-screcching (tha 

latter may also Include speed changes). 

The problaa of daceraining the rpcinir:'! delay 
crajeccorics can be solved by exaaining Che loci for 
fuel consuapcion versus Cine for Che ainiaua fuel 
specified-range (SR) and ainiaua fuel specifled-ciaa 
(ST) probleas plotted in Fig. 8. First observe chat 



Fig. 3 Fual-Ciae loci of optiaua crijaccories. 


Next asaksse that cha spacifiad final ciaa Cj 
is greater chan c^^i. There the fuel consuapcion 
read froa the ST locus ac Cj will be equal to or 
less Chan chut read from the HC locus for all 
d > d^. It will be exactly equal for Che value of 
d ■ d^, naaaly, where Che two loci are tangent co 
asch other. Thus, in order Co niniaiss fuel con- 
suapeion for C;t > Cg,^, the ainiaua distance pach 
of langch d^ should be scratched by the differ- 
ence dj - d^. The detailed shape of Cha path- 
SCretch aanauver is not cricical. However, cums 
should be done at soall bank angles so as co oini- 
aiza bank-angle-induced drag, which was nsglectad 
in this derivacion. 


Tha final fcap in solving ;ha optiaua delay 
problem is co show hou c,j, the aaxi.aua ciaa to 
fly a specified range d|, can be coaputad froa the 
algoricha developed in Lhe preceding section. The 
Haalltonian, Che state equation, cha cnilsa cose c 
and cha transversality condition for tha flxed-tlaa 
problca with specific energy’ as indapandant 
variable are 


H * ain 


j c^Wg - t(E^) ^ - C(Eg) 


''up'”upl 

'^dn’"dn 


iE>o 


|E| 


E<c 


(33) 


d(t 


up 


+ c,„)/dE 




t (E.) 
up i 


'dn'"f> 


Cruise cost: 


C(E^) - Mn CjWj(V^.E^) 5 C,W,,^^(E^) 


Constraint: T ■ P, L ■ W 


(34) 


9 



wh«r« Cyp Is s elM Inscsnc during cllsb, a 
clM inscsnc in descsnc assaursd poslClvc in cha 
backward cIm dlraceloni and eer^p la cha claa 
apant In crulsa. To ylald cha ssm optlBua erajac- 
cory. cha Haallttoniana of Eqa. <16) and <33) mac 
bacowa idantlca,l< This raquiraa Chac ^ ■ -A In 
Cq. (16) ba zaru. or froa Eq, (22), 

C-M. -f C 

\(E.) - win S. , 0 (36) 

® Ve 

AC cha yalua of ■ ^copc ■Inlaicas Cha 

oparand In Eq, (36)# wa auac hava 

Cg - -CfVf (37) 

Buc froa^Bq. (34), Wf > ^^faln' charafora, 

Cc ■ ‘^fWfHin(Ec) * axpraaaion for la 

aubselcutad tnCo Eq> (36), cha ainlalsaclon oparaclon 
doas Indaad ylald A - 0 ac Cha ninlalzing valua of 
Vc. Thus, ac aach cmlaa aiiargy a Claa coac 
coapucad from 

C^(E^) - -aln CfW, (38) 

'^c 

will yield cha naxiaua delay crajacCory, The Crana- 
varsaliey condiCion (Eq. (35)) cranaforaa co cha 
following expreaalon: 

(H + dc(d6jai„<Ec)/«c>) - 0 (39) 

Exaaplea of opclaua delay crsjacCorles coa- 
pucad by cha procedure described hare as wall aa an 
alCamacaclva darlvaclon can ba found In Ref. 8. 

Exaaplaa of Opclavi Tralaccorlaa 

This sacClon prasanca savaral axaaplaa of fuel 
and coaC opciaua CraJacCorias. Addlclonal axaaplaa. 
Including cha affacla of winda and of air Crafflc 
concrol conaCralnca, are given In Raf. 7. 

Hlnlaua Fuel Traiaceorlaa . Figure 9a shows 
exaaplea of alnlaua fuel craJaccorlas for ranges of 
100, 200, and 1,000 n. al. Tha life, drag, and 
propulsion aodals used In Chase axaaplas ara rapre- 
sancaclva of cha Boeing 727-100 equipped wich 
JT 8D-7A anginas. The cakaoff waighc la 150,000 lb, 
winds ara zero, and che acaoapharlc aodal la cha 
1962 ICAO SCandard. For Che 200-n. ml. range, boch 
cha consCralnad-chruse (solid line) and cha frea- 
Chruac (dashed line) CraJaccorlas ara shown. Also, 
for cha 200-n. ml. range. Fig. 9b shows Cha corre- 
sponding alclcuda versus airspeed profiles. The 
conacralned-chrusc crajaccorles for che 100- and 
200-n. ml. ranges concain shore cruise segmancs 
below Cha opclmum crulsa alclcuda of 32,000 fc. 
Opclnua crulsa alclcude Is reached for rangaa 
greacar chan abouc 230 n. ml. For che 1,000-n. ml. 
flighc, che opclmum crulsa alclcuda Increases ac a 
race of abouC 2.5 fc/n. mi., because of fuel bum 
off. Differences bacwaen Che conaCralned- and fraa- 
chrusc craJacCorlaa are apparenC only near che Cop 
of cha climb, where che free-ChrusC CraJecCorlas do 
noc concain a cruise segmanc. The difference In 
fuel consumpclon bccween che conacralncd- and 
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(b) Alrspaed-alclcude profiles, 200-n. mi. range. 

Fig. 9 Minimum-fuel crejeecorles. 

free-chrusc crajeccorles for che 200-n, ml. range 
is 23 lb ouc of a cocal of 4,800 lb. This rala- 
clvely small dlffarance would seam co JuaClfy Cha 
compucaclonally simpler conseralnad-chrusc mods, 
especially In an on-board ImplemantaCion. The alml- 
laricy of cha climb and dascanc profiles for dlf- 
fsranc ranges also offers opporcunlCias for 
simplifying che on-board algorichm. 

Evaluation In Plloced Simulator . Frequently 
the question Is asked how much fuel use and costs 
can be raducad by using on-board fllghtpath opci- 
mlzatlon In airline operatlona. For several 
reasons, chls question Is not easy co answer. 

First, chare doas noc exist a standard refaranee 
trajectory for comparison. Second, there can ba 
significant variations In flight caehnlquas becwsan 
different pilots. Third, It Is difficult to achieve 
rapeatabillty In CraJacCorias because of unknown 
dlscurbancas from weachar and air traffic concrol. 
Reallsclc evaluacions, therefore, should Include 
statistical analyses of savings from s large number 
of trial flights. Hare we excerpt results from a 
limited evaluaclon, using a piloted simulator.* 

The conscralnad-chnisc version of this algo- 
rithm was Implamented on a DC-10 simulator and 
Incagracad with flighc director and autopilot sys- 
tama. Qualified DC-10 airline pilots first flaw 
che slmulacor on a 220-n. ml. fllghC in a manner 
racommanded by chair airline flighc manual. Than 
they rapaaCad the flighc, using Cha flight director 
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CO follow crajoctorlta lonoracod in real ciM by cha 
al)|iorlCha> Thrav rapfaaancaclva crajacci>rlas froM 
chtaa ataulaCor flighta tira pluCcad in FU* 10 ua 
altituda varaua fual cnnauaiad. Nota Chat cakaoff 
And touchdown altiCudaH ora 3,000 fe and 300 fc, 
raapaccivaly, corraapondlng eo cha aleltudaa of 
runwaya at cha particular city palra uaad In cha 
alBNtlacor fllghea. 



Fig. 10 Standard airline procodura coaparad with 
fual- and coat-opcimal crajaccorlaa, 220-n. ai. 
range • 


The fual cpnMuaad and flight tiaa for Cha stan- 
dard airllna procadura ware 9,300 lb and 37 win, 
raapaccivaly. Cruiaa alcituda waa 24,000 fc. Tha 
fual-opcinua trajectory cliabad acsaply and nora 
slowly to 32,000 ft, whereupon it began Inaadiataly 
an Idal-ehruac deacenc. Tha fuel consuaad and cha 
flight CiM wara 8,750 lb and 42 ain, reap^iccivaly. 
The fual saving of 730 lb, or BZ of total fuel, la 
highly significant in airllna econoaiiea whara sav- 
ings of only tZ ara conaidarad laportant. On the 
ucher hand, cha Incraased flight Cine may ba unde- 
ainsblc, but reflects cha face Chat the cost of cIm 
la naaunad to ba zero in cha fual opeiaaim ease. To 
eliminate the time-fuel trada-off from cha conpari- 
aona, a coat-optimum crajactory waa flown with tha 
cime-cost factor aelcctad ao as to achieve cha 
sane flight cIm, 37 min, as the standard sirlira 
procadura. This trajectory required 9,150 lb of 
fuel. Thus, even at fixed arrival cine cha optimum 
trajeccory reducad fuel consunpclun 350 lb, or 3.7Z, 
relative to cha airline standard procedure. Clearly, 
Che fual consumption difference at fixed arrival 
tlma provides a useful measure of cha afflciancy of 
.in airline flight procedure. Fuel savings on 
longer-rangc flights were aimilar in magr.^cuda but 
lass when expressed as n parcenc.tue of total fuel 
consumed. 

Pilots Judged the optimum crajeccorics flown 
with Che aid of a flight director no more difficult 
to fly than the standard procadura. They considered 
the coat-optimization feature an essential element 
of a future on-bourd syscam. Variations of this 
algorithm have been i;icorporaccd in performance- 
msnagaaient systems being built by several avionics 
Ruinufaccurers. 


Tarmlnal-Arua Flithtoath Manaaameni 

Taninal-arsa ( llghcpath MnagaammC problosui 
are concaptuslly and analytically mors difficult 
cibin cha an routa prohloma dlscusaad in cha pracad- 
ind saccion. Thair complaxity artaas from cha 
dynamically mora complex modelK required tor synthe- 
sizing afficianc and flyabla crajaecorlns Tha 
synchaain problem can ba stacud as spacifving an 
algorithm for ganarating a tralectorv from an 
arbitrary initial statu vector (Xi,Yi,hi,Hi,Vj) to 
a final stats vector (Xf,Y(,hf .Hf,Vf • This 
so-callad capture trajectory problem is illust rated 
in Fig. li. Note that in tha I igura cha final point 
lisa on a backward excsnsion ul cha cunway caitCar- 
Una and that the final heading is equal to cha 
nuiway heading. The distance * row the final point 
to touchdown la piloc-spacified Tha capture algo- 
rithm must synchesizs efficient trajectorias 
rapidly and without failure for a wide range of 
initial and final state vectorb 




Fig. 1 1 Tha capture trajectory problsm in cha 
cccvinal area. 

At the present it is not feasible to inplsMnC 
the optiMl control solution for this fivc-scata- 
vsriable problem, as we wara able to do for the 
one-state-variable en route case. Nevertheless, 
optiMl-ooncrol theory played i crucial role In 
deriving the r.lgorithm. It was first used to derive 
the structure end characteristics of sxtrsmum tra- 
JsetorisS (chose that satisfy necasaarv conditions 
of optiMlicy: see Ref. 4) for '^wo reduccd-urdcr 
■ubprobiama of the original (ive-atatc-variable 
problem. Then som of the derived characteriacics 
ware incorporated in the deelgn of the algorithm. 

thus, the first step In deriving a practical 
algorithm is to separate Che synthesis into two 
sssenCially indspendenc problsM chat can be solved 
In sequence.; The first coneisce of synthesizing a 
horizontal plana (two-dimensional) trajectory chat 
connects the Initial position and heading 
to the final position and heading (Xf,Yf,H{). The 
second consists of synthesizing airspssd and alti- 
tuds profiles chat connect Initial and final speeds 
and altitudes (Vj^.h^) and (Vf,hf), rsapactlvely, 
along the known horizontal path. This cschniqus 
was dsscribed In Ref. 10 in connection with a study 
of four dlMnsioiisl guidance and baa been refined 
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••varal Although «oeh aubproblM la 

IndopanJancly optiuliadi coiiblnlng tha trajactoriaa 
froB liha two aubproblaaa doaa not ganarally yiald 
optlMBi aolutiona to tha original probloa* Navar- 
thalaaa, trajactoriaa for tha Boat fraquantly 
ancountarad Initial eondltlona hava baan found to 
glva parfocBanea raaaonably eloaa to tha optlBUB. 
Moat laportantly, pllotad alaulatlona and flight 
taata hava ahown that tha afflclaney of tha coBputad 
trajactoriaa axeaada that of pllot-ganaratad flight 
patha. 


takaoff and cllab-out profllaa. Claarly* the opcl- 
BlaaCion of auch changaa In aircraft configuration 
la an Intagral part of profile aynthaala. 

Both optlBlaatlon of tha configuration and tha 
aynthaala of afflclant profllaa can ba hand lad with 
a raforaulatad varaion of tiin anargyatata nodal 
davalopad for tha an routa caae. 

Tha anargy rata, Cq. (6), can bo wrlttan In 
tha fom 


Svnthaala of Horlaontal Fllaht Patha 

Horlaontal plana trajactoriaa ara clioaon to 
alnlBlaa tha length of the path or, aqulvalantly, 
the tlna to fly at conatant alrapaad and aaro wlndn. 
Thla problan can ba foraulatad aa a nlnlnuB tlaa, 
optlnal control problaa with atata variables X, Y, 
and H, and with bank angla ^ as tha control 
variable. Analysis shows chat the axtranua trajac- 
corlaa conalst of an Initial turn, followed by 
althor a straight line or another Cum In tha oppo- 
alca direction, and a final tum.^' Turns ara flown 
at naxinuB bank angla (^nax)* alnlnuB Cum 

radius Rgiin la conputad fron the relation 

*.ln " ♦.ax 

where Vg la tha ground apaad. A repraaantatlon 
of the solution In teraa of nanauvarlng regions 
covering tha entire chraa-dlBenalonal atata apace 
la given in Ref. 14. 

Tha coBputaclonally Boat efflclant solution and 
one chat has baan laplamantad In a flight ayataa la 
baaed on a set of closed-fom equations for coBput- 
Ing all possible extraanai crajectorlas. For each 
aat of Initial and final conditions, the algorltha 
coaputas up to six different extraanaa trajectories 
and chan chooses the one with the ahortaat path 
length. The alnlatn Cum radius for each cum la 
specified saparataly In the algorltha, using 
Eq. (40). Since airspeeds and ground apaada will 
not actually raaaln constant, an estiaate of the 
luxlBUB ground speed la used In Eq. (40) to ensure 
chat the chosen Cum radius will not cause the bank 
angle Halt Co be exceeded. A conaarvatlva eaciaate 
for Che aaxlmua ground spaed la the algabralc aua of 
Che aaxloaBt alrapaad In Che Cura and Che wind speed. 
Note chat astlaaClng tha aaxlaua lapllas a weak 
Incaracclon between Che assuaed Independence of Che 
horizontal and vertical synchesls probleaa. 


« _ dh . V ^ 

dt 3T i dt 


(41) 


By using tha relation dh/dt a Vy In Eq. (41) and 
dividing Eqa. (6) and (41) by V, wa obtain two 
axpraaslona for a quantity E„, defined aa tha 
nonaalized energy race 


■ (T(r)cos(a + v) - 0(a,6f,V))/W (42) 


1 dV 




(43) 


with L x W as a constraint. Tlia cos(o * v) fac- 
tor accounta for tha possibility of chruat vector- 
ing. Equations (42) and (43) provide a siapllfylng 
dlchotoBy In tha profile synthesis. Ac a particu- 
lar alrapaad the (noraallzed) anargy rata la first 
datamlnad froa Eq. (42) by choice of appropriate 
control's. Including thrust T, Chruat-vaetor angla 
V, flap angla 5f, and angle of attack a. Then , 
the choaan anargy rate bacoaaa an Input Co Eq. (43), 
which decaralnes tha alrspead-alcltuda profile by 
apaclfylng either y or dV/dt. 


Since a convancional aircraft has three con- 
trols and a typical V/STOL aircraft hsa four to 
aehlava a specif lad energy rate, wa hava an axeass 
of one and two controls, raspactlvaly, over the 
alnlama nuabar naeded for a unique slmltanaous 
solution Co Eq. (42) and L " V at a given altltuda 
and alrapaad. These extra degrees of fraedoa In 
the controls can be exploited to nlnialze thrust 
and, charafore, fuel flow at every specified energy 
rate, alcltuda, and airspeed. This optlnlzaclon 
problaa la restated In equivalent fora as the aaxl- 
aizatlon of energy rata for a given thrust settings 


i (ir) • laax 


Tha horizontal capture algorltha, coabinad with 
contlnuoua display of the trajectory on an elaccronlc 
nap, has received aany favorable coiBMncs by fllght- 
cest pilots. Furthamora, it can be used as a sub- 
prograa In algorlthas for coaputlng variable-radius 
turns and crajeccoriaa through a aaquanca of way- 
points. For Cheoa reasons Che derivation of Che 
algorltha la Included In the appendix. 

Synthesis of Alrapaed-Altltuda Profiles 

The doalnanc feature of the alrapeed-alclcude 
profiles lA Che temlnal area la the deceleration 
and descant segaanc eo achieve landing approach 
conditions. In this aegaant, cha flaps are extended, 
and, for certain types of V/STOL aircraft, the chruat 
is vectored to Increase drag and to coopansate for 
dlnlnlahed aarodynaadc life at low speed. SlaUlar 
changes, except In reverse order, occur in Che 


wlch constraint L(n,v,a,Sf) > W. The aoxlalzatlon 
Bust alao obey various inaquallty constraints on tha 
controls, auch aa Halts on the flaps, the angla of 
attack, and thrust-vector angle. Wa nay Interpret 
the aaxlalzaclon operation as a tachnlqus for col- 
lapsing cha Bulclple controls Into Che single 
vsrlabla 

Tha lower bound of energy rata attainable by 
choosing tha controls by Che process of Eq. (44) la 
the BaxlBua energy rate atcalnable at olnlaiBi- or 
idle-chnisc setting. (In a V/STOL aircraft, at low 
speed, tha alnlnuB chruat may be higher chan Idle.) 
Hmevar, aore negative anargy rates. If the aircraft 
can attain chon, oust also ba generated by the 
BOdel in order co encoapaso the entire flight enve- 
lope of Che 'aircraft. For in aircraft without 
thrust vectoring, such rates are obtained slaply by 
Increasing the flap angle b^.yond the opclaua 
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obt 4 ln«d from Eqt (44), For 4 V/STOL mlrcraft with 
Chrxiat vaccorlnn uha problM Itt mora complaK, In 
chla caaa, cha Lntruduucion of a aacond critarlun, 
ahich adnimitaa ttiu thruat daflacclun anRia for 
aacli apaclflad anorgy rata, halpa to ganarata cha 
raaalndar of tht .iccaiiubla anargy rata raglon.'’ 

In a practical Imptamancatlon of chla procadura 
for a vactorad-chruac STOL aircraft, cha functional 
ralationa baevaan cha anargy rata and cha apctmiiad 
cuncroLa ara pracompucod and acorad in cha on-board 
t'anputav in ■ultidimanalonal Cablaa.^* Intarpoia* 
Cion of ancriaa in Cha cablaa dacarminaa cha opClaua 
concroln (or apaclflad vnluaa of anargy rata, air- 
apaad, altltuda, timbianc camparncura, aircraft 
waighc, and bank angla. by nondlaanaionalliing and 
combining aoma of thaaa variabLaa, cha acoraga 
raquiraamnc for cha cablna waa roducad to only 
60U worda, avait for cho ralaclvaly complax vacCorad- 
Chruae STOL aircraft, Tha acoraga raqulramant would 
ba conaidarably loaa for u convanclunal Jat crana- 
porc altcrafc. A almpliflad varaion of Chaaa Cablaa 
for Cha STOL alrcriifc uxampla la givan in Fig. , 12 
aa ploca of maxlaaiin and minianim anargy raCaa, 
and ^nmln v«raua cha oqulvalanc airapaad, 
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Fig. 12 Envalopa of aitccgy raca (or a vaccorad- 
chruac STOL aircrafCl W - 38.000 lb, aaa lavai, 
■candard Camparatura. 

After davaloping tha aachod for aalacclng Cha 
uuncrola, wa can now prucaad wich Cha aynchaaia of 
the profilaa. In principle , chia probltmi ia idanci- 
cal to cha fixad-ninga, an routn problam pravioualy 
Holved. Howavar, cha shore rangaa and variacy of 
opar.-»clonal conacr.iincx characcarlatic of canainal- 
iiraa flight Juacify 4 ximpllfiad approach baaad on 
nntching cha ganarnl characcoriatlca of optimum fuel 
airapaad-alcicuda profilaa. Wi< briefly explain Cha 
racionale for chla machod with rafaranca to daacanc, 
which la tha moat difficulc caar. 

Hinimum fuel daacanc crajaccoriaa auch aa Choaa 
llluacracad in Fig. 9 ara characcariiad by a mono- 
conic dacraaaa in anargy from cruiaa until tha 
apaclflad final (landing) airapaad apd alcicuda ara 
nac/ Furchamore, at carainal-iiraa alcicudaa below 
10,000 fc, cha crajaccoriaa firac daacand at 
approximaCaly conacanc indlcacad airapaad of 
240 knoca (for a boaing 727) and idle chroccla to 
Cha landing alcituda. Tha daacanc ia (ollowad by 
rapid deca location Co cha landing apaad ac nearly 


conacanc alcicuda. Howavar, ubacructiona and noiaa 
raaericCiona along Cha approach path usually do not 
panaic laval-fllghc dacalaracion ac low alcicuda. 
Inaeaad, cha dacalMr.iCion ia parlurmad in ahallow 
daacanc or in lavai I light at altlcudaa batwaan 
1,000 and 2,000 ft .ihova cha runway. 

To provida (laxibility in cha ahaplng of cha 
profilaa during aimuUanaoua dacalaracion and 
daacanc. familiaa of dacraaaing (and by axtanaion, 
incraaiing) energy profilaa are ganaraCad aa a 
function of two panimacara, o and c, Tha firac 
paramacar, o, salact.s tha fraction of ■inimum 
(Mxiaum) avAilabla anargy rata, Knaln* <^Ban^ 
bo uaad (or daeraaning (Incraaalng) anargy. The 
valuaa of "nux pracomputad and 

acnrad at each indicaca~air apaad. aa pravioualy 
axplalnad. Tha aacond paramatar, c, datanainaa 
Cha fraction of the xalacCad anargy rata Co be uaad 
(or dacalaracion (accalaraCion) Than, for parcicu- 
lar chuicaa of a and c, Cha anargy race, airapaad 
rata, (llghcpach angla, alcicuda raca. and ground 
apaad ara computed >it) followa: 

K “ -''^nmin • " ^ <*5> 

V - g« E„ . 0 S 1 i 1 (46) 

Y - U - c)E^, (47) 

h • (48) 

A • V coa Y (49) 

Note Chat ^nmin ^ ^ * providaa dacalar- 

acion vichouC daacanc ; >. ■ 0 u'ovidaa daacanc wich- 
ouc dacalaracion: and 0 < c « providaa wlaul- 
canaoua dacalaracion and daacanc. Hinimum fuel 
daacanc ac alCiCudaa nearly Che aama aa landing 
aicltudaa uaually requires following the 
contour in cha anargy-rnca cablaa. Thua, cha opci- 
mum value of o Is unity. Howovar, for aoma air- 
craft, .such aa cha vactorad-cht uat STOL type man- 
Clonad previously, Cha optimum value may yield 
anargy racas coo nagaclve (or safe operation. A 
value laaa chan 1 ia also naca^iaary to raaarva a 
mar^.in for c losad-loop control along cha computed 
path. A prttccical upper limlc for o ia about 0.9. 
Furchaneoi'e, maximum dacalaratlun and daacant llmlca 
are also enforced during profile synchaals. 

Tha Btructura of cha profl lea ia modeled after 
chat in Fig. 2, axeape tliac an additional conacralnt 
la impoaad. An aircraft flying In Cho Carminal area 
la ganarally noc allowed co climb above ica initial 
approach alcicuda for tha purpoaa of opciaiiBlng 

cha trajaccocy. Rjichar, it muat hold this altitude 
uncil starting cha daacanc Cor landing. Howavar, 
while flying at alcicuda h^, it may change to a 
fuel-af flclenc carminal-araa airapaad V^, For cho 
STOL aircraft, V(- ^ 140 and for a Jac cranupore, 
auch aa a 727 3 240 knoca IAS. Another cri- 

terion (or chooaing is co maat a apaclflad 
landing clma, aa required in four-dlmanaional 
guidance applicaciona. 

Tha variouB rulaa above can now ba cpeibinad 
co ganaraca complete profilaa. The synchaals bagina 
with tha backward time intagraclon of Eqa. (46), 
(46), and (49) from final condiciona hf. Vf, uaing 
cha apaclflad 0 and c. If tha alcicuda raachaa ica 
cargat value of h(^ before Cha airapaad raachaa lea 
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target value ot we aet c ■ 1 and then con- 
tinue the backward tlaa Integration until the air- 
apead haa alao achieved Ite target value. Whan 
netting c the fllghtpath angle la forced to 
lero and the energy rata la uaed entirely for accal- 
aratlng (in backward tlaa) toward V^. On the othar 
hand If the airepead reachaa Ita target value before 
Che alcltuda does, we aac c ■ 0. Thla scope the 
airepead change and cauaea the energy rate to be 
uaed entirely for Incraaaing the altitude coward lea 
target value of h^. Whan the accond and last 
variable reachaa Ita target value, wa set a •• 0, 
Chat la. fen • 0, thus conpletlng tha backward tlM 
integration. Next, wa begin a forward cine Intagra- 
tlon troa «-.ha currant aircraft position to gat tha 
distance re<iulrad to change epead tram to 
with c ■ I. Let the dletancas for the backward 
and forward Integrations bo dj,t) and d^f. raspoc- 
tlvaly. K valid trajectory haa been ganaracad If 
the cruise dlacance d^. coaputad froa 
d(. - df - dib - dif. Is nonnegatlva, that la, If 
dc < 0. If dc la nsgatlve, the synthesis has 
failed, because the aircraft Is coo close to the 
capture point P(. For purposes of on-board lapla- 
aentatlon, tha laportant feature of this algorlcha 
la that It synthaalses an efficient trajectory in a 
single Integration. 

Figure 13 lllustratea the various ssgaents of 
a synthaaised approach trajectory for a STOL air- 
craft, with a ■ I and c - O.S. We aaauae for 
nlapllclty chat Ennln " -0>13, a conatant. Other 
paraaetera are indicated In the figure. Note chat 
the Initial descent at y ■ -7.5* reduces to 
Y ■ -3.7S* In order to allow the aircraft to dacal- 
arata to tha landing apeed of 100 ft /sac. 



Fig. 13 Exaaiple airspeed-altitude profile for a 
vectorad-chniac STOL aircraft. 


Tha algorlcha alao corrects tha alrspaad dacal- 
aration for known wind ohaars, which are coaputad 
froa a knowledge of Vy(h) If available. The wind- 
shear correction factor la 

- -(dV,,/dh)VY 

which la added Co Che right side of Iq. (46) Co 
obtain the wlnd-ahaar-corractad alrapaad rate. 
Furtheraora, the raferanca controls are corrected 
for tha efface of nonsero bank angla on Induced 
drag by aodlfylng the weight In Che lift equation 
aa follows: L ■ W/coa 4> 

Xncagratlon tlaa-ntaps vary during aynchaals. 
During turns, dacaleratlons, and acealaraclona It 
la 1 aac; during altitude changes at fixed spaed 
and heading It la 5 sac. Total tlaa for synthesis- 
ing a cosiplace crajactory consisting of a conplex 
horlioncal path, such aa that shown In Fig. 16 In 
the Appendix, ai^ an alrapaad-altltuda profile, 
alnllar to the one shown In Fig. 13, is about 2 sac 
on tha particular airborne conputar used In a 
recant flight axperlsMnc. This coaputer haa an add 
tlaa of 6 Msec and a stultlply/dlvlde tlaa of 
24 usee. Whan the trajectory synehasla is claa- 
sharad with navigation and other nacasaary coaputa- 
Clona, the coaputlng tlaa Incraasaa to about 6 sac. 

Noncircular Capture Traiectorlaa 

Tha coaputaclonal slapllclty of the preceding 
algorlcha depended on the a priori aaparetlon of 
the ayncheaia problaa into (nearly) Independent 
solutions for Che horlioncal paths and alrapead- 
alcicuda profiles and on Che choice of slaple gao- 
aecrlc foraa for horisontal paths. Racently. 
Kralndlar and Nauaan^^ atudled fual-optlaun capture 
trajactorlaa under leas restrictive conditions. 

For crajacCorles containing a fairly long straight- 
line aegaenc between the Initial and final tutna, 
they found that Che axtreaua Crajectoriaa had 
approxlaacely the following characteristics. As 
before, the airspeed In the stralghc-llne segaanc 
tends to be close to the silnlxua-fual-per-unic- 
dlatance speed. However, the final Cum la flown 
at ■axlsan bank angla while the aircraft Is decel- 
erating to the landing apeed. Equation (40) shows 
Chat Che resulting turn la a spiral of deersasing 
cum radius. For Curas with large heading changes 
('ISO*) this saves fuel by reducing the tlaa and 
distance flown. 

The horisontal capture algoriChai for circular 
turns can be used In aavacal s:;aps to generate 
approxlsuitely conscant-bank-angle spiral turns, as 
shown In Ref. IS. The nachod la illustrated In 
Fig. 14. For slapllclty, asauas chat the initial 
cum has fixed cum radius Using Eq. (40) 

calculate at Pf for a bank angle a few 
degrees 44 lass chan 4^^ ^g ~ ^f ^wi> 
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Fig. 14 Illustrating conacrucclon of noncircular 
final cum. 



wh«r« !■ the wlntl-apaad coaponant In tha final 

haadi>rig diracelon. Than calculaca cha circular 
horiaontai path with final cum radtua Rf ■ R^. 
ivViiC tncagraca cha apaad, altlcuda, and ranga aqua- 
Cionai Eqa. (46) . (47). (40) • and (49) backward 
froa ?( until tha ground apaad la auch Chat 
4 • 4aax radlua R^. Thla ceaplacaa tha flrat 

atap. At chat point compuca a naw and larger turn 
radlua Rj ualng V. > V] 4- whara Vj and V,f| 
ara cha airapaad and along-crack coaponanta of cha 
wlnd-apaad vactor at tha and of the flrat atap. 

Than apply tha clrcular-cum algoricha a aacond 
clna and raauaa cha backward incagracion. Conclnua 
atapplng tha cum radlua until alchar cha crulaa 
apaad la achiavad or tha baginning of tha final turn 
ia raachad during tha backward Incagracion. In cha 
axanpla of Fig, 14, four auch atapa wara nacaaaary. 

Obvloualy thla horiaontal path algoritha 
raquiraa aonawhat nora conputatlona for aach ayn- 
chaaia. Tha apaad of cha availabla on-board con- 
pucar will datamina If it can ba uaad in an 
application. 

Svatan Inolawantation and Tliaht-Taat Raaulta 

Tha on-board inplaaancatlon of tha ayacaai ia 
baat'd on two nodaa of oparatlon. In cha flrat noda, 
rafarrad to aa Cha “predicclva*' ooda, naw trajac- 
coriaa ara aynchaaiiad ona aftar cha ochar aa 
rapidly aa poaalbla. Upon conplatlon of aach ayn- 
chaaia, Cha ayacan chacka Co dacanalna whachar cha 
piloc haa callad for cha aacond, or "Crack" noda Co 
ba angagad. If auch ia cha caaa, cha pradlcClva 
noda ia carttinatad. Than tha noac racanCly ayntha- 
aiaad craJacCory ii raganaratad and crackad in caal 
clna by a cloaad-loop guldanca law. 

Tha piloc activataa cha pradlctlva noda by 
aalecClng a waypoint to bo capCurad on a flxod cra- 
JacCory, which la praacorad. Synchaala of cha cap- 
cura crajactory baglna aftar tha navigation ayaton 
haa conpucad Cha currant poaitlon and valocicy con- 
ponanCa of cha aircraft. Tha flrat atap in Cha 
synchaala conputaa cha horizontal Crajactory paran- 
i^cars; ic ia always succaaaful. Tha aacond scop 
conputaa cha airapaad and alclcuda proflla; sona- 
tinas ic can fail. For axanpla, if cha horizontal 
path ia vary short and cha dlffaronco in anargy 
baewaan initial and final poaltlona ia largo, a 
flyabla trajectory along chat path nay not axle. 
Howavar. a failure to synchoalzo ia unlikaly in land- 
ing approaches initiacad aavaral nllaa fron cha 
capcura point, Che usual siCuaclon. If Ic occurs, a 
diagnoaClc naasaga ia dlaplayad Co Cha pilot indi- 
cating Cha reason for cha failure. Aftar a failure, 
Che algorlchn auconaclcally rapaata cha aynthoals 
proCesa, ualng updacad position and valocicy vectors. 
Also, cha piloc can always correct tha failure to 
synthesize by flying tha aircraft away fron cha cap- 
ture point or by selacclng a nora distant waypoint. 

Synchaalsad crajaccorloa ara dlaplayad to tUa 
piloc on a oMp-llke cachoda ray display callad a 
nulcifunccion display (MFD). Figure 15 gives an 
exanplo of crajaccprlaa dlaplayad on Cha MFD. Tha 
piloc haa salacced waypoint 3 on Cha fixed Crajac- 
tory (dram solid) aa tha capture waypoint. Tha 
daahad-Una crajactory starting at aircraft poslclon 
Pj indicates to tha pilot that a valid capture 
trajectory has bavtn synchealted. In cha exanpla tha 
pilot did not engage tha crack nods at P^ buc 
instead flew tha aircraft in Che direction of Pj. 
Beewean Pj and P. the dashed capcura crajactory 



Fig. 15 Horizontal flightpacha dlaplayad on nulcl- 
f unction display. 

was rafraahad approxinataly avary 6 aac. With tha 
naw ganaraclon of airbomo conputara, Cha rafrash 
rata, which is dataralnad by cha apaad of synthasls, 
can ba incraasad co a nora desirable onca-par-sacond 
race. Whan cha aircraft raachad P} cha piloc 
alaccad Co engage Cha track noda, causing cha last 
capcura Crajactory co ba frozen and radram aa a 
solid line on cha MFD. Ac chat tlM. cloaad-loop 
Cracking of tha frozen capcura crajactory was 
initiated. 

It should be r.ocad chat in Cha Crack noda, Cha 
aynchaslzsd crajactory ia not rafraahad, chough 
this nay ba desirable if winds or cranslencs in 
navigation introduce largo cracking errors. The 
aofewara of cha on-board operating syatan haa bean 
configured co add chin capability in cha future. 

Closod-loop cracking la parfornsd by a parcur- 
baclon guidance law spaclflcally daalgnod co oparaca 
In concert with Che output of the crajactory ayntha- 
ala algorlchn. Parcurbacion stacea in the feedback 
law include errors in airapaad, altitude, flightpath 
angle, and cross-crack position. Controls aro 
chruac, thrust angle, pitch, and roll angle. Tha 
feedback law was deolgned with cha help of quadratic 
opclnal synchaala cachniquas and thus differs in 
aavaral raapaeta from conventional aucopllocs. 
Datalla of cha doalgn for a vaccored-chruac STOL 
aircraft are given in Raf. 12. 

A powarad-lift, vac corad-ch rust STOL aircraft, 
rafarrad to by NASA aa cha Aupnentor Wing Jet STOL 
Raaaarch Aircraft, was salacCad as the Cast vehicle 
for evaluating the covseept.'^ This aircraft Is 
aqulppsd with a general purpose digital conpuCar 
and flexible navigation and control ayacazM. The 
complex aarodynaalc and operational characcerlacics 
of this aircraft preaancad major challengaa in cha 
design aa wall aa opporCunltlaa for damonstraclng 
cha value of automacad flightpath managamanc. Ona 
cast aaquance compared cha fuel conaiaipclon of a 
aynchaalzad and a pilot-flown approach crajactory. 
Both approaches began at Cha same initial dlatance- 
co-couchdown (40,000 fc), airapaad (140 knots), and 
altitude (3,000 fc). Tha plloc-flown approach was 
made with the aid of a flight diraecor syatan. Tha 
aucomaclc approach consumed 301 lb of fual, and cha 
manually flown ona consumed 500 Ib. The difference 
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bcewian chaa ta conaldarad fairly rapraaancatlva 
ovar a wlda ranga of initial condiciuna. 


Concludina Kaaarka 

Trajactory optialaation occura aa an aaaantial 
atap in tha daaign of on-board flightpaeh aanagoBant 
ayataaa. Optiaal control theory providad tha nacaa- 
aary analytical fraaawork for gaining inaight into 
tha charaecariatica of officiant aolutlona and for 
unifying divaraa problana. Hovavar. tha darivatioo 
of practical on-*board algorithaw vaa nora atrongly 
dapandant on phyaical raaaoning rind on Judieioua 
ainpllficationa of thti problaa than on axaet Inpla- 
nantation of optimal-control thaory. In practica. 
tha critical toot of an algorithm ia not whathar it 
ia optimum in a mathamatlcally praciaa aanaa, but 
whathar it can conaiatancly outperform ita compati- 
tor, who in thia eaoa ia tha unaidad pilot. Simula- 
tion and flight taata hava ohown that both tha 
an routa and tha tormina 1-aroa algorithm moot thia 
criterion. 

Navarthalaaa, olnco tha tochni^uaa daaeribad 
hara ara firat-ganaratlon aolutlona » opportunitiao 
ara abundant for further Improvananta in parfommnea 
and for automation of other difficult pilot taaka. 
Chief among such taaka ara following behind an air- 
craft at a apoclfiad minimum dlatanno and marglng 
smoothly into a otrsarn of aircraft on landing 
approach. Attention also naada to bo given to 
Integrating tha an routa and taminal-araa flight- 
path aanagamant algorithms. 

Tha potential for computar-dlractad trajactory 
aanagamant in military applicationa la widely racog- 
noiod, but on-board algorithms chat outperform an 
experienced pilot ini for examploi typical air- 
coabat situaciona are far aore difficult to obtain. 

In Che near tarn, pronlelng areas for applying the 
approach developed hare ara in autoaatod guidance 
of renotaly piloted vahiclaa and in noncombr.t flight- 
path aanagamant. 


Appendix 

Tha expreaslona for synches ixing horiaontal 
capture crajactorlee for flying an aircraft from a 
given Initial position and heading to a epacifiad 
final poaicion and heading in a minimum diatanca 
are derived below. 


Tii im-Straiaht-Tum Tralactoriea 

'The turns ara arcs of clrclas and tha straight, 
portion of tha trajectory muse ba tangent to tha 
Initial and final circles. Since tha Initial and 
final turns nay be althar clockwlae or counterclock- 
wise, Chare ara four possible combinations of turn- 
ing dlractionsi two with the initial and final turns 
in the Sana dlractlon and two with Cham in oppoaica 
directions. Figure 16 llluscrataa one solution of 
each type. If a given pair of clrclas is aneiraly 
separate. Chat is, if no part of one circle lias 
within tha other, it la possible to draw four tan- 
gent lines batwaen the pair. However, vector 5 
along the tangent line from the Initial to tha final 
circle coincides with the direction of rotation at 
both tangent points for only one of tha four tangent 
lines as shown in the figure. 



Fig. 16 Turn-scralght-curn cases. 


In fig. 16, the final poaicion and the origin 
of tha coordinate syacam ara located on tha runway 
cancarllna; however, the derivation is Cor arbitrary 
locations. Furchansora, all variablaa ara dafinad 
so that tha darlvacion appliaa to all posaibla 
combinations of turning dlractione. 

figure 16a ia for Che caaa in which both tuns 
ara in tha same dlractlon, and the tangent vector 
5 does not cross Q; in fig. I6b, Che cuma are in 
oppoaita directions, and 6 erosaaa Initially 
the aircraft is at (Xi,Y^> in soma inertial Car- 
tesian coordinate aystam %rith heading dafinad 
as poaiciva clocbriaa from the X-axia, and 7x 
a unit vector in tha direction of tha velocity. 

The vector distance from (Xj^.Y^) to the canter of 
the Cum la given by OxR} where <*** radius 

of cum, and vector normal to 

and positive to Che right of Thamfora, the 

vector from (Xj^tY^) Co the cantar (XCx,YCi) ia 
RjOx ■ right cum and -R^Ox for a left turn. 
The directions of cha cuma are accounted for by 
writing the radius vector as RxSxOx, whara $x " '^1 
for right turns and Sx “ -1 for laft turns. Simi- 
larly, cha direction of tha final cum Is denoted by 
Si. 

Tha aircraft aovaa along tha circle from 
(Xi,Yi> CO the tangent point (Xj,Yx), which l«a a 
radius vector RxSxUi. Tha tangent vector from 
(X, ,Yx) at tha and of cha initial cum to (X|,Y,) 
at cha baginning of cha final turn ia D. The 
radius vector at (Xx,Yx) is RiSjU], but since 0] 
and Ox must ba normal to 6, U; • 0,. Likewise, 
cha headings Hj and H, at Cha two tangent points 
arc equal. Tha final turn ends at (Xf,Yf) with 
heading Hf and radius vector RjSxO^. 

Using this notation wa can writa 
D + RjU^Sx “ RjUjSx ^ Q 
or 

q • 6 + UjCRxS,, - RjSx) . (Al) 
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•ndt charaforat alnca 5 and 0, nra parpandlcular • 

D s |BI • /Q* - ^l,Si - 

whara by daflnltlon 

Q * |(}| - /<XC, - JtCjVii ♦ <YC, - YCii* <A3) 

(t can ba aaan frun Iq. (A2) that no raal aolutlon 
axlaca If Q < |R;S{ - RiSxl • Whan tha tuna ara tn 
oppoalca diracctona. Si ■ -S|« and thara ia no raal 
aolutlon for Q < di Rx)f that l 0 i If tha clrclaa 
Ihtaraact. On tha othar hand> for retatlona in tha 
tana dlraetlona, Sj • and a raal aolutlon axlaca 
unlaaa Q < |lt} - Ri|> that la. unlaaa ona elrcla 
llaa anelraly vlthln tha othar. Fron laonacrlc 
conatruccion It can ba abown that thara alwaya axlat 
at laaac two raal aolutlona. Fron tha daflnltlon of 
tha radlua vactora, ona can wrlta for tha raal 
aolutlona: 


RjOiSj 


/-RjSj am 

\ R^Sx coa M^/ 


and 


/xci ■ h\ 

i“iSi - ( ) 

\xci - <1/ 


Equating Eqa. (A4) and (AS) glvaa 


Similarly, 


XC| ■ Xx - R(Sj aln H^, 
YC, - Y, RjSx cot 


XCj " Xf - RtSx aln 


YC 


2 * ^ f *" Xjb 2 coa Hj. 


(AS) 


(A6) 


(A7) 


Tha radlua vactora at tha tangant polnca can ba uaad 
In cha aaaw nannar to conputa tha conponanta of 
S; and X| 


X, « XC, ^ R^Sj ain Hj 
Yj - YCj - RjSj coa Hj 


»,02 


(A8a) 

(A8b) 

(A9a) 

(A9b) 


D coa H, - (XC. - XC,) * (R.S, - R.S,)aln H, 
D aln K. - (YC, • YC,) - (R.S, - R,S,)coa H, 


(A12) 

IquaClona (A12) can ba aolvad (or cha tangant of 

H,: 

(YC, - YC,)D - (R,8, • R,S,)(XC, - XC,) 

*'* ■ (kd, ~ xd,)D (n^Sj - R,S,)fW, - YC,) 

(A13) 

Iquaclona (Ad)-(A9> and Iq. (A13) eonplataly apaclfy 
a eapcura crajactovy for any conblnatlon of S, and 
Sj. Howavar, the langth of cha trajactory la alao 
naadad In ordar to dataraina which of cha faaalbla 
Crajactorloa glvaa cha nlnlaun dlatanca, Tha flrat 
cun angla la 


(AA) whara 


TR, - (H, - Hj) 4- 2wC,S. 


{ 0 if S,(H, - Hj) c 0 
1 If S,(H, ■ Hj) . (I 


(AU) 


{ 0 If 

I If S,( 


(AIS) 


and cha aacond turn angla la 

TR, - (Hj - Hj) 2wC,S 

whara 

,(Hf - H,) i 0 
5,(8, - H,) 0 

Finally, cha total langth of cha eapcura path la 

dj - |6| -fRjlTRj * RjiTR, (Al«) 

Tun-Tun-Tun Tra1 actor laa 

Thara ara at noat four Cun -Cun-Cun- typo 
Crajactorlaa for the horltontal-captura problan; 
nanaly, two aach of Cha rlght-laft-rlght and laft- 
flghc-lafc pactana. Howavar, wo allnlnata cwo of 
than by raqulrlng chat cha nlddla cun axcaad n 
radlana . ^ * 

Tha problan la llluacratad In Fig. 17 for Cha 
rlghc-lofc-rlghc pateon. Tha vactora In cha figure 


Subtracting (A8a) from (A9a) and (A8b) fron (A9b) 
glvaa tha componenca of 6: 


X, - X, • XCj - XC, + (RjSj - R,S,)aln H, 


Y, - Yj - YCj - YC, - (RjSj - R,S,)coa H, 
Another axpraaalon for Cha conponanta of 5 la: 


(AlO) 


X, - X, - 0 coa H, 

Y, - V, • 0 ain H, 


(All) 


Equating Cha corraapondlng palm In Eqa. (AlO) and 
(All) glvaa 
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•rt d«(ln«d as before , reeognlilng that in this case 
$2 ■ S|. To satisfy ths requireaent that the Middle 
turn exceed r radians, its canter aust be on the 
opposite side of 0 froa the straight segMnt of 
the righC'Straight-right solution shown for coaparf* 
son, Furtharaore. no thrae-arc solution axists for 
q > R, -f 2R,. 

Let Hq be the heading of Q, defined in the 
text, and define a unit vector with heading 
angle, H, as follows: 

»# - Hq ♦ *1 I 

Than V, is perpendicular to end points in the 
dlraction of flight idiera Q intersects the clrela 
of the initial turn. Fora the law of coainaa 

Q* ■¥ (R. + R,)* - (R. ♦ R,)* 
coe A' - 2Q (R^ R,} 


where 


where 


Tl, 


c. - 

TR, 


Ct ■ 


> (H, 

- Hg) 

+ 2wC,S^ 

0 if 

(H, - 

Hg)8 t 0 

1 if 

(H, - 

H,)S < 0 

■ (H, 

- Hg) 

- 2xCtSx 

fO if 

(H, - 

H,>Si i 0 

[l if 

(H, - 

H,)S, < 0 


(A27) 


(A28) 


The total length of the trajectory is therefore 

dj • RjTRj + RjItR,, + RjITR,: (A29> 


coe B' 


Q* + (Rj + R,)* - (Rj + R,)* 
2Q ^R, + R,) 


(AIB) 


Tha direction of turn is accounted for by defining 


The length of trajectories for all feasible 
pairs of $2 and Sj is coaputed and tha trajectory 
with the shortest length is selected. A FORTRAN 
listing of the algoritha and sons applications are 
given in Ref. IS. 


A - SjA' 


<A19) 


References 


and 


(A20) 


Using these definitions it can be seen froa the 
Fig. 17 tfiac 


Hj - H, + A 1 
H, - H, - B + ir j 


(A2l> 


Froa the definition of H, and the equations derived 
previously 


sin H, ■ 


cos H. ■> 


SiCTCj - TfCj) 


-SiCXCj - XCj) 


(A22) 


(A23) 


Equations (A21) can be used with double-angle trig- 
onoaetrlc identities to obtain expressions for 
sin Hj, cos Hj. sin H, , and cos H«, in teraa of Hg, 
A, and B. Changing appropriate subscripts in 
Eqs. (A8) and (A9) gives 


, - XC 2 -*• R^Si sin H,1 
] - YCj - R^Sj cos Hg| 

, - XCj + RjSj sin H, \ 

2 — ^^2 " ^2^1 COS Hg I 


(A2A) 


(A25) 


The turn angles are calculated as follows 

TR 2 - (J»5 - Hj) + 2ffC,S2 (A26) 


where 


0 if (H, - Hj)S^ i 0 

1 if (H, - H,)S2 < 0 
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